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I. Introduction 

Recently the mode number of the Dirac operator on the lattice has been used to determine the 
mass anomalous dimension for both confining and conformal theories [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 

II, 12, 13]. This method has been shown to work well in the case of confining theories, where the 
running of the mass anomalous dimension can be followed from the perturbative regime through 
to the confining regime by fitting to different ranges of eigenvalues and combining results from 
several lattice spacings [8]. For QCD this has recently been done with the addition of a continuum 
extrapolation, which gives very good agreement with perturbation theory [11]. For theories with 
an IRFP the method gives very precise values, see e.g. Ref. [6], although the dependence on the fit 
range and bare coupling is less clear. As stated in Ref. [10] it would clearly be desirable to apply 
the method to a theory where analytic results are available for comparison. 

In this work we apply the mode number method to the rif = 2 massive Schwinger model. The 
Schwinger model has a long history as a toy model for lattice simulations, see e.g [14, 15, 16, 17, 
18, 19] and references therein. This is a very attractive toy model for testing the mode number 
method, as in addition to having an analytic solution for the mass anomalous dimension, it lives 
in only two dimensions and so is inexpensive to simulate numerically. This was recently done 
in Ref. [10], but only a single gauge coupling was considered and the expected value for y* was 
not found. Here we show that by measuring y as a function of the eigenvalue and using several 
different bare couplings, as advocated in Ref. [8], it is possible to follow the running of the mass 
anomalous dimension over a range of energy scales from the UV to the IR, where we find a value 
which approaches the analytic result. We also apply the method to the rif = 0 theory which has a 
non-zero chiral condensate, where we find the expected behaviour in the IR. 


2. Mode Number Method 


In a mass-deformed conformal field theory (mCFT) in d dimensions, the chiral condensate 
goes to zero with the mass as 

(2.1) 

where y, is the mass anomalous dimension at the IRFR 

The spectral density p(co) of the Dirac operator at small eigenvalues ft) also goes to zero with 
the same exponent [1,2] 

d _i 

lim lim p(ft>) ft) 1 '<■ . (2.2) 

m—>0V— 

In addition, for an asymptotically free theory in the UV limit this equation can be written in 
the same form but with y* replaced by y(g|), the 1-loop mass anomalous dimension which is a 
function of the renormalised coupling g\ [8], 


d 


lim lim p(ft)) ft) 177 * 8 * 1 
m —>0 V— 


-1 


(2.3) 


The idea is then to fit the measured spectral density to an equation of this form to extract an 
effective mass anomalous dimension y(A) as a function of the eigenvalue A. This quantity is equal 
to y* in the IR limit A —> 0, and to the perturbative 1-loop y(g\) in the UV limit. 
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In a theory with a non-zero chiral condensate, such as QCD, we expect instead p (ft)) = p (0) / 
0 in the chiral limit. Note that if we nonetheless perform a fit to the form of Eq. (2.2), this would 
give “y*” = d — 1 in the limit X —> 0. 


3. Massive Schwinger Model 


The massive Schwinger model consists of n/ Dirac fermions of mass m which live in 1+1 
dimensions, and interact through a U(l) gauge field. The euclidean Lagrangian is given by 

= ¥f( x ) + 8^n (*)) + m] \f/f(x) + -F^ v (x)F^ v (x), (3.1) 

where the index / runs from 1 to rif. For rif > I massless fermions, the spectrum of the theory 
consists of one massive boson of mass p = sjn.fg 2 /n , and rif— 1 massless bosons [20]. 

For rif < 1 the theory has a non—zero chiral condensate, whilst for nf > 2 the condensate 
vanishes in the chiral limit. Assuming the fermion mass is light compared to the U(l) boson mass 
(m < p), and the volume is large compared to the U(l) boson mass (pL 1), then the mass- 
dependence of this condensate is known in two limiting cases [21, 22]. 




m for mL <C 1 

m (n f -\)/{n f +\) for mL ^ 1 


(3.2) 


Comparing this to Eq. (2.1), we see a very similar picture to the one we have for mCFTs; in 
the UV limit we have a free theory with y = 0, whilst in the IR limit y —> y, = I /rif. Hence for the 
rif = 2 case we consider here, y* =0.5. 


4. Simulation Details 

The implementation on the lattice used here follows Ref. [14], using the plaquette gauge action 
and unimproved Wilson-Dirac fermions. Configurations are generated using the Hybrid Monte- 
carlo (HMC) algorithm, with trajectory length z = 1, and the number of integration steps tuned to 
keep the acceptance rate in the range 60 — 90%. A(2000) configurations were generated for both 
rif = 2 and iif = 0, at f5 = 5,2,1 on lattices of size 32 2 . A smaller number were generated at stronger 
bare coupling, /3 = 0.5,0.1. Each configuration is separated by 20 HMC updates, and on each the 
lowest 300 eigenvalues are calculated using Chebyshev accelerated subspace iteration [23]. The 
Wilson-Dirac operator explicitly breaks chiral symmetry, so at each value of /3, fc = 1/ (2m + 4) is 
tuned to K ~ K c (determined here as the value which minimises the lowest eigenvalue aQ. o). Some 
additional configurations were generated with larger masses, and on smaller lattices, to check for 
finite mass and finite volume effects. 


5. Fit Function 

On the lattice we measure eigenvalues O 2 of the massive hermitian Dirac operator M = m 2 — 
lj)~ . These are related to co as co = yO 2 — m 2 . In general the lowest eigenvalues will be affected 
by finite volume and/or finite mass effects, and the highest eigenvalues will be affected by lattice 
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p(H) for n f =2 



n 


Figure 1: Log-log plot of eigenvalue density p (Q) Qi -r versus Q for n j- = 2 and various values of /j. 
The slope changes from y ~ 0 at large eigenvalues (red line shows the analytic prediction of the slope in the 
UV limit to guide the eye) to y ~ 7 * = 0.5 at small eigenvalues (black line shows the analytic prediction of 
the slope in the IR limit). This demonstrates the importance of the fit range in determining 7 *. 

artefacts, but we expect Eq. (2.2) to apply for eigenvalues in some intermediate range (Ojr < ft) < 
cOuv■ The strategy is then to determine y(ft>) for many values of ft) within this range. We consider 
two ways to fit the data. The first is to fit the measured mode number directly, which assuming 
finite mass and finite volume effects are negligible can be done with the simple fit funtion 

v(aQ) = A ( aQ ) l +r( aSl ) , (5.1) 

where A and y(aQ) are free fit parameters. Alternatively one can first construct the spectral density 
by binning the measured mode number, then do a fit to 

p(aQ) =B(aQ)m^~ l , (5.2) 

where B and y(aQ) are free fit parameters. 

In both cases we determine y(aQ) by fitting the data in the range aQ. ± A, and both fits should 
give the same value for 7 * at the IRFP. Using Eq. (5.1) gives more precise values for the fitted 
parameters because no binning of the data is required, but the price to pay is that finite volume 
effects below the fit range still affect the fitted parameters. Using Eq. (5.2) avoids these finite 
volume effects, but increases the statistical errors, and there may also be a systematic dependence 
on the bin size. 

To investigate a wide range of energy scales we use several different values of /3 = 1 /g 2 , which 
correspond to different values of the lattice spacing a. The resulting y(apQ) can then be combined 
by rescaling the eigenvalues aQ in terms of the lattice spacing at some fixed value of /J [ 8 ], in this 
case we use ap =5 0 . In the Schwinger model a °< 1/y^/3, so we define the rescaled eigenvalues as 

/ nr\ l+ y( a p a p'> 

Q = a p QpU^\ . (5.3) 
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7(H): n f =2, Mode Number Fit y(H): n f =2, Spectral Density Fit 




Figure 2: Mass anomalous dimension y vs eigenvalue for the nj = 2 theory. Left: mode number fit, right: 
spectral density fit. Different colours correspond to different values of /3 . The results appear to be consistent 
with the analytic prediction of y —► 0 in the UV (large eigenvalues), and y —► y* = 0.5 in the IR (small 
eigenvalues). 


6. Results rif = 2 

Fig. 1 shows the spectral density p(Q) as a function of Cl on a log-log plot. The slope of the 
curve determines the value of y, which runs from y~ 0 at large eigenvalues to y ~ y, = 0.5 at low 
eigenvalues. For each value of f5 the points have been multiplicatively scaled to make them lie on 
a single curve, this is only to make the plot clearer to the eye, and does not affect the slope of the 
curve. This plot shows the importance of considering a range of scales when determining y using 
this method. 

Fig. 2 shows our determination of y(Q) for the rif = 2 theory. The left plot shows the values 
determined from fits of the mode number to Eq. (5.1), while in the right plot the spectral density, 
constructed by binning the mode number data, is fitted to Eq. (5.2). Several bin sizes are used in 
the plot to check that there is no systematic dependence within errors on the choice of bin size. For 
each value of y, the x error bar shows the range of eigenvalues used in the fit, and the y-error bar 
shows the statistical error of the fitted value. 

At large eigenvalues, both fits give values of y consistent with zero, the analytic prediction in 
the UV limit. For intermediate eigenvalues the two fits give different values, but as one goes to 
smaller eigenvalues the two start to converge, as they should, and appear to be consistent with the 
analytic IR prediction of y, = 0.5 at zero eigenvalue. Going to smaller eigenvalues would require 
a larger physical volume, either by using a larger lattice, or by increasing the lattice spacing. Note 
that although the fit ranges have been chosen such that finite volume and finite mass effects are 
negligible, no continuum limit has been taken so lattice artefacts remain. 

7. Results rif = 0 

We also apply the same method to the quenched theory with nf = 0. Fig. 3 shows the resulting 
determination of y for this theory. Since this theory has a non-zero chiral condensate, we expect 
the mode number method to find “y*”~> d — 1 = 1 in the IR, which is consistent with what we find. 
The rif = 0 and rif = 2 theories clearly have different values of yfor small eigenvalues, however at 
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7(H): n f =0, Mode Number Fit y(H): n f =0, Spectral Density Fit 




Figure 3: Mass anomalous dimension y vs eigenvalue for the nj = 0 theory. Left: mode number fit, right: 
spectral density fit. Different colours correspond to different values of . In the IR limit (small eigenvalues), 
y = 1 corresponds to the expected chiral symmetry breaking. 

intermediate and larger eigenvalues there is little qualitative difference between these two theories, 
which shows the importance of considering a range of energy scales when trying to determine y, 
for theories which may or may not have an IRFP. 

8. Conclusions 

For the rif = 2 massive Schwinger model, we were able to follow the running of the effec¬ 
tive mass anomalous dimension over a wide range of energy scales using the Dirac mode number 
method, from y ~ 0 in the UV to y ~ y, in the IR. Our determination of y at the lowest measured 
energy scale is within ~ 10% of the analytically known value in the IR limit. This accuracy could 
be improved by going to lower energy scales, either by using larger lattice volumes, or by using 
stronger values of the bare coupling to increase the lattice spacing, and performing a continuum 
extrapolation. A weakness of this method is that to avoid finite mass and finite volume effects, one 
is always forced to work at non-zero eigenvalue, and the appropriate form of the extrapolation to 
zero eigenvalue is not known, so that this distance from the IR limit remains as a systematic error 
that is difficult to quantify (unless one already knows the answer as is the case here). 

We also applied the same method to the rif = 0 case, which has a non-zero chiral condensate. 
Whilst we found the expected behaviour at small eigenvalues, for intermediate eigenvalues the two 
theories have very similar values of y, despite their vastly different IR behaviours, which underlines 
the importance of considering a range of scales when determining the mass anomalous dimension 
using this method. 

Flaving confirmed that the Dirac mode number method reproduces the analytic predictions 
for this toy model, it would be interesting to use it to investigate the extrapolation to the IR (zero 
eigenvalue) limit, as well as the extrapolation to the continuum limit. 
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